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Abstract 

Using a multipole expansion of the radiated field generated by a classical electric 
current, we present a way to interprete the bremsstrahlung spectra of low energy 
heavy ion collisions. We perform the calculation explicitely for the system 
at 84AMeV and compare the result with the experimental data of E. Grosse et at 
Using simple model assumptions for the electromagnetic source current we are able 
to describe the measured data in terms of coherent photon emission. In this context, 
the information contained in the measured data is discussed. 


1 Introduction 

The angular distribution of bremsstrahlung emitted from accelerated charges strongly 
depends upon the structure of the source current. Correlations between the dynamical 
properties of the charges show up in the spectrum of radiated photons. A multipole 
expansion of the radiated held |]^, |^, ^ has the advantage, that special structures like 

symmetries of the current enter directly in the coefficients of the expansion. Dominant 
structures of the current give dominant terms in the expansion, structures and shapes of 
the spectrum provide therefore information about the current. At least one can conclude 
from structures in the spectrum, that the current is somehow structured as well. Using 
the multipole expansion and simple model assumptions for the electromagnetic source 
current as well as for the geometry and the dynamics of a heavy ion collision, it is the aim 
of this paper to show that the experimental data of E. Grosse et al. can be understood 
as the coherent radiation from a correlated motion of the nucleons. In the following we 
show how the expansion is obtained. In section we discuss the current that will be used 
for the calculations. With this current the exact and the approximated spectrum are 
calculated and compared in section ||. Finally, in section § further assumptions, needed 
to describe a heavy ion collision, are made and included in the calculation, which is 
explicitely performed for the system at 84AMeV. 


1 


Starting from Maxwell’s equations, assuming a periodic time dependence of the sources, 
and looking for a solution of Kirchhoff’s differential equation in spherical coordinates one 
finds a representation of the electric and the magnetic field in terms of Bessel functions 
and vector spherical harmonics. Due to the duality of E and B in Maxwell’s equations 
one gets two different solutions for E and B whose linear combinations are hnally the 
general expressions for the helds. 

More generally one can expand the vector potential, i.e. the wave function of the 
photon, in two transversal and one longitudinal vector field, spanning a 3-dimensional 
space at each point. The vector fields are designed in a way that they represent outgoing 
waves in the far held. Since a real photon has only two polarizations one only needs the 
two transversal helds in the free space outside the sources. 

Thus, the magnetic held reads |]^: 

B(r,u,) = kY. (1) 

Im 


{k = \ k I),whereas in the far held one gets for the electric held 

E{f,u!) = B X e(r2) . 

The functions 0^^^’‘\k,r) are dehned by 
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with Oi the outgoing spherical Bessel function. The coefficients of the expansion of B in 
terms of the read @] : 


a 


VA{i) 

m 


(cn) 


dvr 

T 




(4) 


and obey the relation 
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The (k, f) are dehned in the same way as the above with the regular spherical 

Bessel function 

k'T 

Ei{kr) = —^ j e)d cose , (6) 

Pi{cosd) being the Legendre polynomials. The expressions for the coefficients can be 
simplihed by rewriting the current j{r,uj) in terms of its Fourier transform and using 
Rayleigh’s plane wave expansion for the exponential function 
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Throughout the calculation we used the orthogonality relations 
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The complex conjugate of reads j^*{u}) = (— 
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2 The current 


The classical 4-current of a charged particle moving on a trajectory x(r) in proper time 
can be written as 0 


f{x, t) = - x{t)] \t=xO{r)= e J dT^5‘^[x - x(r)] . (10) 

Rewriting the (5-function by its Fourier integral one finds, that the current for a colliding 
particle can be understood as the Fourier transform of 


ji^(k) = -ie 
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assuming that the collision can be described by a simple kink in the trajectory. Coulomb 
deflection effects, which are important in low-energy collisions, are neglected here. Gen¬ 
eralizing o to the case of several particles suffering several collisions we obtain 


fik) = -ieY.Y. 
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where the indices count vertices (i) and particles (j), respectively. Xij in the equation 
above accounts for the space-time history of the scenario. 


In the following we will make a specihc model for the current, considering four colliding 
particles with equal charge and mass. These particles undergo two subsequent and inde¬ 
pendent two-body collisions. The hrst pair interacts at the space time point = (0, (0)) 
and the remaining two collide at = ( 2 / 0 ,^)- For simplicity both collisions are assumed 
to share the same reaction plane and absolute value of the deflection angle. The Fourier 
transform of the current in the center of mass frame then reads: 
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with 

r(«. 6) = , - 

1 — ax -\- by 1 — x^ cos ifk 

a = j3\\, b = and x = cos 6k. For sake of simplicity we only consider a temporal distance 
At = yQ between the second collision and the hrst collision. The current ( p!3D has therefore 
been constructed in a symmetric manner. For At = 0 it is not only invariant under a 
parity transformation but also mirror symmetric with respect to axes both parallel and 
transverse to the beam direction. These symmetries are broken with hnite At. 


3 The spectrum 


The spectrum of bremsstrahlung photons radiated by charged particles is given by the 
simple expression H 
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with dk = The angular dependence of the spectrum is determined by the scalar 

products in the denominators of (0) = poko — \p\\k\e{flp) ■ e{flk), e{flp) ■ e(Qk) = 

cos 6p cos 6k + sin 6p sin 6k cos((pk — <Pp) • 
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The spectrum of radiated energy of the process ([T3| ) therefore reads (/ = ujN) 
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Inserting the current (jl^) and averaging over the azimuthal angle (pfc, leads to 
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with the abbreviation Q{a,b,x) = y^( 1 + ax)2 — b‘^{l — x2). 
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The photon spectrum is given by the Poynting vector, i.e. 
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After changing the order of summation Z) = S S Z Z and averaging over the 

Im l=0m=-l m=—oo 1=0 

azimuthal angle (p, i.e. averaging over all possible reaction planes, terms containing the 
product of two with different m drop out, which yields: 

dl _ _ 

duidilk 27r Yn 

We will describe the spectrum with the nonvanishing terms of lowest order in Z. In a 
similar way the photon emission was calculated for nuclear collisions below the Coulomb 
barrier 0. 




( 18 ) 


The occuring products of the vector fields determine the angular depen¬ 

dence of the spectrum. Due to the averaging over ip these products are real and depend 
only upon cos Ok- 

We use the far held limit for the functions O: 
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Both, the product and respectively, yield the well known angular 

functions P| 
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For equal I the interference terms are asymmetric functions of cos 9k- The evaluation 
yields for I = 2 
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and for I = 1 (only m = ±1) 
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For symmetric charge distributions, i.e. currents with even parity, the coefficients of these 
terms vanish. 


The interference of dipole {I = 1) and quadrupole (/ = 2) contributes with two different 
kinds of terms: both, the interference of magnetic and electric components of B and 
the interference of magnetic and magnetic as well as electric and electric components, 
respectively. The latter is asymmetric in cos 9k @ 

= d);)(2)*(|fc|,r) ■d;^W(|fc|,r) 
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Hence, for a current with even parity the interference reduces to the product 
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Figure 1: Comparision of the exact bremsstrahlung spectrum and the spectrum calculated 
by the hrst terms of the multipole expansion for a scattering angle 9p = 3.6°, the initial 
rapidity y = 0.2 and the hnal velocity f3f = 0.5/9^, for At we chose the large value of 
20fm/c. The spectra are symmetric with respect to 9k = tt/2. The transverse scattering 
generates hnite radiation aX 9k = /2 even for soft photons (/ = 2, m = 2 moment), the 
temporal phase between the two pairs of interaction causes the oscillation of the different 
terms in the expansion, i.e. the increase of the / = 2, m = 1 moment and therefore photon 
radiation in beam direction. 


Here and in eq.(^) we used that we restricted the motion of the particles to a plane. 
The scalar product of with in eq.(|^) becomes a scalar product with eq, or Cy. 

Choosing, without loss of generality, the x — z-plane as scattering plane, the coefficients 
aQ^''\u}) vanish. For symmetric collisions the coefficients and vanish, 

too. The possibility to substitute by \l\ex further simplihes the calculation since 
now complex conjugation of the coefficients and the switching of the sign of the 

index m decouple (see eq. (|])). 

A comparision of the exact spectrum (pT]) and the hrst terms of the multipole expansion 
(^) in Figure |l| shows the quality of the restriction to terms of lowest order in /. This 
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restriction is valid as long as the initial and final rapidities can be described well enough 
by 

y-0 + Y ( 25 ) 

For higher energies, terms of higher order in I have to be taken into account, i.e. the 
bremsstrahlung spectrum of ultrarelativistic particles is composed of infinitely many terms 
of the multipole expansion and their interference. Although the approximative character 
of this treatment is obvious, the multipole expansion offers a very instructive way to 
understand the behaviour of bremsstrahlung emission. 


4 Comparision with data 


To compare the spectrum obtained with the multipole expansion with experimental data, 
one needs to make further assumptions. We will perform an analysis of the radiation 
measured in the system at 84AMeV by Grosse et al who determied the 

photon emission in minimum bias collisions. One therefore has to average the spectrum 
over all impact parameters b, weighted with the geometric cross section of the collision^: 


da 1 /■2’r fbmax 

- = — dip / bdb - 

dudO, u Jo Jo dudVt 


(26) 


bmax can be chosen as 2R + x, where R is the radius of the colliding nuclei and x is the 
surface thickness of the mass distribution. 

We treat noncentral heavy ion collisions in the way proposed in [Q, i.e. the charge 
Z{b) involved in the collision is obtained from the geometric overlap integral P] 


V{b) = J d^xe{R^ -x'^-y^- z^)e{R^ -y‘^-{z- bf) (27) 
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V ■ 


V and Z are the volume and the charge of the nucleus, respec- 


and hence Z(b) = Z- 
tively. 

Earlier calculations exhibit the dipole structure^ of the photon spectrum, assuming ei¬ 
ther the validity of incoherent summation due to an uncorrelated motion of the nucleons 
or contributions from nucleus nucleus collisions at large impact parameters whose 
radiation should possess a dipole-like structure as well Even a mixture of coherent 
and incoherent radiation has been suggested for the description of the radiation [Q. In 


the following we want to propose as an alternative mechanism radiation produced by the 
correlated, collective motion of the nucleons. The current (|T^) is considered as a coherent 
elementary process in the heavy ion collision. For that reason we constructed the cur¬ 
rent (|T^ in a symmetric manner since the averaged final state of a symmetric heavy ion 


^One assumes that, as it is valid for soft photons, the cross section of the radiation process factorizes 
into the cross section for the scattering process times the photon spectrum. 

^According to the remark given in the following of eq.(^^ the spectra of high-energy collisions should 
neither be called dipolar nor quadrupolar, unless these names refer to the first nonvanishing electric 
component of the expansion 










collision is symmetric as well. Eq. (|^ represents the simplest current which serves our 
purposes. We account for all processes which can be described by (JT^) considering one to 


four charged particles. For collisions involving neutral particles the corresponding terms 
in (|l^) are simply set equal to zero. We further assume that all nucleons collide at the 
origin, i.e. the current simply has to be multiplied by a temporal phase and integrated 
over time. The resulting spectrum of the nucleus nucleus collision is then similar to that 
of the nuclei feeling a box-like force ~ 6{T /2 — t)6{t + T/2) in the time interval T, which 
here represents the time of the collision. This time dependence of the collision is also 
obtained, when one treats the collision of the two nuclei in the frame of a shock wave 
model, neglecting the spatial distance of the two outward travelling shock fronts, i.e. one 
performs the calculation as if the two shock fronts were sitting at rest in the origin (with 
this assumption, radiation is calculated in [p!3| ). 

Adopting the shock wave picture, we can estimate the minimal collision time with a shock 
velocity jSgh equal to the speed of light. 


T = ^ 

min ( n , i \ 

7i(A +1) 


(28) 


D is the diameter of the nuclei. For noncentral collisions the time is taken to depend 
linearly on the diameter at 6/2. 

To obtain the correct behaviour of the spectrum with respect to the photon energy, 
one can introduce a phase space correction |T^ 


R2{s) X^R{s,mi,m2) s 


(29) 


with s = s — 2ujy/s and X{x, y, z) = {x — y — zY — iyz the kinematical triangle function; 
mi and m 2 are the masses of the colliding particles, respectively. This way of treating the 
spectrum of high-energy photons is strictly justihed only for the incoherent summation 
of the spectrum. Since no adequate formalism is available at the moment, however, we 
adopt this correction factor (^^ also for coherent emission, i.e. we multiply the resulting 
spectrum by (^^. To describe the high-energy tail of the spectrum it is unavoidable to 
consider either three- and four-nucleon clusters which feed their energy to the photon pro¬ 
duction in a collective manner, as proposed in |T^. We increase the nucleon momentum 
by the mean Fermi-momentum which has to be taken into account since it is of the same 
order of magnitude in the considered system. Consequently, the available phase space of 
the photons gets enlarged. 

In additionally to the decreasing available phase-space, the time dependence of the colli¬ 
sion determines the emission of high energy photons, too. Since the energy distribution 
of the spectrum is the square of the Fourier-transform of the acting force, the spectrum 
of the considered process is poportional to 


dl 

dujdVt 


cx 
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sin^(a;T) 


(30) 
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Since lim = x both, uj and T enter symmetrically in the spectrum, one obtains 

a;—>0 * 

the soft photon limes for either a; —0 or T —>^0 (for all photon energies). Therefore, 
for short collision times one obtains a nearly constant, flat spectrum. In this way, the 
collision time controlls the slope of the spectrum for not too large photon energies. A 
lower bound for T, however, is given by ([^) and has the approximate value of 5fm/c. 

As mentioned earlier, the presence of asymmetric collisions can be expected,which are 
the consequence of fluctuations of the charge distribution in the nuclei. The probability 
for these asymmetries of the charge distribution in the overlapping volume as a function of 
the impact parameter we parametrize in the following way: For nearly central collisions it 
is zero (the probability for a symmetric collision is 1) and smoothly drops to 1/2 when the 
mean charge of one nucleus in the overlapping volume becomes < 1. The probability for 
symmetric collisions approaches 1/4, the remaining fourth is the probability for symmetric 
but radiationless collisions, i.e. neutron neutron collisions (see Figure]^). 



Figure 2: Parametrization of the probability for asymmetric and symmetric collisions as 
a function of the impact parameter. For convenience, the functions used here are for both 
cases a hyperbolic tangent. For symmetric collisions: Psymm = ^|tanh(artanh(a)-((6/i?) — 
1)) + I and for asymmetric collisions: Pasymm = ^x^anh(artanh(a) ■ {(b/R) — 1)) + |. a 
determines the stiffness of the parametrization and is taken to be a = 0.99. 

As possible asymmetric and symmetric collisions we consider all possible combina¬ 
tions of up to four charges in (|^ with one exception: instead of the collision of only two 
charged particles (one from the projectile and the target) possessing a temporal distance 
At between the scattering processes, we consider the symmetric case with two particles 
at each case, since it is unlikely that all particles scatter in only one transverse direction. 
The nonvanishing coefficients are, besides a factor two which is canceled since only half 
of the charge enters in the latter case in the calculation, in both cases the same. In the 
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latter however only coefficients with even m exist. 

The ’’delay”-parameter At breaks possibly present symmetries of the scattering process. 
In this way At controlls (with respect to a;) the strength of the contribution of different 
terms of the expansion. Processes with a symmetric charge distribution in the overlapping 
volume of the colliding nuclei can therefore (for certain photon energies) generate a purely 
dipolar structure in the photon emission. This observation is similar to the result in . 


The asymmetric collisions generate also the interference terms of types (piD -(P^. How¬ 
ever, since the probabilities of the asymmetries are invariant under a parity transforma¬ 
tion, these terms cancel when averaging over all possible asymmetries. 

A quantitative evaluation shows that the influences of different parameters counterbal¬ 
ance to a certain extent. These are e.g. the collision time and the cluster size for higher 
photon energies and At (for intermediate photon energies) or the strength of the stopping 
(for soft photons), respectively, and the ratio of symmetric to asymmetric collisions. 



0 20 40 60 80 100 120 140 160 180 

^CM 


Figure 3: The bold line shows the spectrum obtained with the dipole (/ = 1) and 
quadrupole {I = 2) terms and their interference. Contributions of collisions with an 
assumed asymmetric or symmetric charge distribution in the nuclear overlap region are 
shown as well. The angular distribution of the data 0 (symmetrized with respect to 
6 = 90°, open squares) can be described quite well. The thin line represents the result 

m 


of an intranuclear cascade model 
between 50 and 100 MeV. 


The cross section shown here is integrated over oj 


Figure ^ shows the spectrum obtained with the discussed assumptions (full line). The 
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dashed line represents the spectrnm calcnlated in the frame of the intrannclear cascade 
model [10|. The angnlar distribntion is not of dipole type and represents the interplay of 


the different mnltipole components and their mutnal interference. Confronting the data 
with the calculation, the measurement seems to indicate an even stronger contribution of 
the quadrupole. One may expect that the photon angular distribution gets modihed when 
a larger number of interacting particles is included in the source current to allow for other 
asymmetries with a different ratio of the I = 1 and I = 2 terms. With better knowledge of 
the charge fluctiations one can improve the assumed impact parameter dependent ratio 
of asymmetric and symmetric collisions. E.g., on has to estimate the influence of the 
repulsive Coulomb forces on the protons in the nuclei. 

One possible signal of the slope of the energy spectrum in Figure ^ is, that the time 
of the collision is much shorter than one might infer from the naive estimation, simply 
regarding the time the nuclei need to pass the distance which gives approximately 

25fm/c for the considered process. Figure ^ shows an arithmetically averaged spectrum 
for collision times between 6 and 16fm/c. The gap at soft photons could be closed by 
adjusting the ratio of symmetric to asymmetric collisions. 


^^C(84AMeV)+^^C 



Figure 4: Comparision of the measured energy distribution at 6^ = 90° with the result of 
the multipole expansion. 

According to Fig. ^ and the considerations made above, the angular distribution 
should change with photon energy. To produce a photon with high energy, one needs 
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many particles which deposit their energy in one photon. Since the probability to form 
large clusters depends on the number of involved particles, collisions with small impact 
parameters, i.e. more symmetric collisions, contribute more to the high energy tail of the 
spectrum. 

For the calculation we used the following values for the parameters: For the mean total 
momentum we assumed 320 AMeV. The hnal velocity was taken to be 8% of the initial 
velocity. At is approximately one third of the actual, b dependent collision time. Since 
the hnal velocity is very small, the effects of the transverse scattering are small as well. 
The scattering angle was choosen to be 14°. These values gave the observed agreement 
with the data. 


5 Conclusion 


The multipole expansion of electromagnetic radiation as already presented in Q] applied 
to bremsstrahlung provides a transparent understanding of the radiated helds. Compared 
to the full calculation the isolated handling of the multipole components and their mu¬ 
tual interference allows for a deeper insight in the inhuence of different parameters and 
properties of the current as, e.g., transverse scattering and the presence of symmetries. 

The fair agreement of the calculation presented in this paper with the data is no 
proof for the correctness of the model. It demonstrates, however, that a completely 
coherent treatment of bremsstrahlung is able to describe the data. The ’’exponential” 
shape of the energy spectrum, which was interpreted as being caused by incoherent photon 
emission from the uncorrelated motion of particles in a thermalized nucleon gas m can 
alternatively be described by the coherent photon emission of the strongly correlated 
motion during the stopping phase in the collision. Our simple model assumptions have to 
be supported by an extended calculation, where a detailed stopping mechanism, as e.g. in 
II, is assumed and the spatial extension of the colliding particles is taken into account. 

We conclude, that the measured data, especially the integrated cross section, do nei¬ 
ther show information about thermalization, though the angular distribution was even 
interpreted as representing isotropic radiation 0, nor provide a proof for incoherent 
photon production. 
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